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and Precession in Accretion Disks 
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0\ : ABSTRACT 



^ . A geometrically thin, optically thick, warped accretion disk with a central source 

of luminosity is subject to non-axisymmetric forces due to radiation pressure; the 



resulting torque acts to modify the warp. In a recent paper, Pringle (1996) used a local 
analysis to show that initially planar accretion disks are unstable to warping driven 
^ by radiation torque. Here we extend this work with a global analysis of the stable 

' and unstable modes. We confirm Pringle's conclusion that thin centrally-illuminated 

accretion disks are generically unstable to warping via this mechanism; we discuss 
■ the time-evolution and likely steady-state of such systems and show specifically that 

. this mechanism can explain the warping of the disk of water masers in NGC 4258 

, and the 164-day precession period of the accretion disk in SS 433. Radiation-driven 

fH I warping and precession provides a robust mechanism for producing warped, precessing 

Qh| accretion disks in active galactic nuclei and X-ray binary systems. 

o 

. Subject headings: accretion disks - instabilities - galaxies: individual (NGC 4258) 



stars: individual (SS 433) 



1. Introduction 

An optically thick accretion disk which is illuminated by the central source of luminosity 
will reradiate the intercepted flux essentially normal to the local disk surface. As first noted 
by Petterson (1977), if the accretion disk is warped {i.e., has a non-constant inclination and/or 
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orientation of the line of nodes as a function of radius) , then the resulting radiation pressure force 
will be non-axisymmetric, leading to torquing of the disk and modification of the warp. 

The origin of this radiation torque can be seen as follows. Consider an annulus of an accretion 
disk with a central source of radiation. If we warp this annulus slightly, it will be illuminated 
by the central source, and, if the disk is optically thick, it will absorb the incident flux. The 
momentum associated with this absorbed incident flux is purely radial, and so the force due to 
the incident radiation will produce zero torque. If the disk is optically thick to re-emission of 
the absorbed radiation, however, the net re-emitted flux will be normal to the local plane of the 
disk. The radiation pressure acting on the irradiated surface due to the re-emitted radiation will 
therefore exert a torque at each exposed point of the annulus. If the annulus is warped, i.e., the 
local angle of tilt varies with azimuth, then the net torque on the ring will be non-zero, because 
the deposition of radiation momentum per unit area of the ring will be non-axisymmetric. In 
general, the net torque will be non-zero provided that there is a radial gradient in either the tilt 
angle /5 (in which case the ring will precess) or the angle of the line of nodes (in which case the 
torque will cause the tilt to change). Note that the disk must be optically thick both to absorption 
of the incident radiation and to re-radiation of the absorbed flux. If the disk is optically thin to 
absorption, then the absorption rate will be essentially constant through the annulus, and the flux 
re-radiated from the two faces of the annulus (the irradiated face and the opposite face, which 
does not see the central source directly) will be the same, so that the re-emitted radiation will 
exert zero torque. If the disk is optically thin to re-emission, then the re-emitted radiation will be 
isotropic, and again the torque will be zero. 
I 1 

Petterson (1977) calculated the influence of this radiation torquing on the disk shape. 



assuming an initially warped disk. However, Pringle (1996) has shown that even initially planar 



disks are unstable to warping due to radiation pressure, so that warping should be generic in 
centrally-illuminated accretion disks. 

Pringle (1996)' (hereafter P96) studied the disk dynamics using a local (WKB-type) analysis, 
in which the variation of the radiation torque parameter and viscosity with radius R was ignored. 
However, such an approach is only marginally justified, since the unstable modes in Pringle's 
analysis have radial wavenumber k ~ In this paper we extend the work of P96 in a global 

analysis; we show that this is in fact analytically soluble, with both stable "normal" modes and 
growing modes. In section § we solve the linearized equation for the evolution of the disk tilt, 
while in section ^ we compare our results to P96 and discuss the time evolution and likely steady 
state of such a system. We show in section ^ that this instability provides a plausible explanation 
for the warp observed in the Keplerian disk of water masers in the nucleus of NGC 4258, and for 
the 164-day precession period in SS 433. 
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2. Analysis of the Disk Dynamics 

Following P96, we define (3 to be the local angle of tilt of the disk with respect to the 
normal to the equatorial plane, while 7 describes the orientation of the line of nodes with respect 
to a fixed axis in the equatorial plane. The equation for the evolution of the local tilt vector, 

t) = (sin /3 cos 7, sin /3 sin 7, cos /3) was derived by Pringle (1992)| (his equation 2.5). We assume 



an isothermal, gas pressure-supported disk sX R ^ Rg {Rg is the Schwarzschild radius), for 
which the surface density S oc R~^^^ and the radial velocity = z/iO'/O, where primes denote 
derivatives with respect to R; vi is the usual disk viscosity, while we denote the viscosity which 
acts to reduce the disk tilt by z^2- This simplifies the evolution equation considerably (c/. P96), 
and, assuming that /? ^ 1 and [3' R ^ 1, it can be written compactly as an equation for W = Pe^'^ , 

dW dW d /I dW\ 

^riT7T + lT75 0^^21777 (1) 



dt dR dR\2 "OR 

where the radiation torque term is 

r = - (2) 

12'kT.R'^Q.c ^ ' 

with L the luminosity of the central source. We make the assumption that the disk is isothermal, 
i.e., it has the same temperature at all radii. Then, if the usual Shakura-Sunyaev viscosity 
parameter a is constant, we find that the viscosity vi scales as z/i = Vo{R/ RoY^'^ where Rq is an 
arbitrary fiducial radius. We further assume that the ratio r/ = vijvx is constant. Since the surface 
density S = M/Svri/i, we can express F as 

r = — Yp- Yp^R = ^oR (3) 
2V2RV' Rl'^ 

where we have assumed the luminosity is powered by accretion; the radiative efficiency e = LjMc? . 
Equation (1) can then be written as 



dt " dR 2 jfj'^ dR V dR 



2.1. Particular Solutions 

For illustration we consider two particular solutions of equation (4). 

a) The steady-state solution, in which the disk has constant structure in an inertial frame, 
can be obtained by setting dW/dt = 0. The general solution for W{R) is then given by 

W{R) = c, + C2 r R''"^ exp ( dR , (5) 
Jo \ r]Uo ] 

where c\ and C2 are arbitrary (complex) constants. Defining 



z 



R'/' = t^{R/RsY'\ (6) 
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(note that z is pure imaginary) this can be written as 



W{z) = ci + 4 



oo pis 



hi/ — as 

Z Jz S 



(7) 



where c'2 is arbitrary. Thus the steady disk has two independent solutions. The first, W =constant, 
corresponds to the trivial solution of a flat disk misaligned with the axial plane. The second 
solution is (formally) singular at the origin (W ~ R^^l"^ as ^ 0). More importantly, we 
require that the viscous torque exerted by each annulus of the disk on its neighbor, given by 
G oc E^/^dW / dR^ is finite as r — > 0. The second solution therefore requires a source of angular 
momentum at the origin. 

b) The 'constant-torque' solution is given by 

W{R,t) = e^^°'l'^R-^l'^ . (8) 

This is the solution that results when the viscous torque G is independent of radius. In this case 
the disk as a whole precesses with angular velocity To/2. The nodes all lie in the same radial 
direction, and the disk is stationary in a frame which rotates with angular velocity To/2. This 
comes about because the radial dependence of W insures that the precession rate induced at each 
radius is independent of radius (c/. equation [3.10] in P96). 

This is actually a special case of a more general class of solutions, W{R,t) = g^^ot/aj^-i/a ^ 
where a 7^ 0. The difference is that for these solutions rj = v^/^x is not constant, as vi oc 



2.2. More General Solutions 

In order to investigate more general time dependence, we Fourier transform with respect to 
time and write dW /dt = iaW; a in general has both real and imaginary parts. We again make 
the substitution (6) to rewrite equation (4) in the form 



This is Kummer's equation {e.g., Abramowitz &: Stegun 1964| , (AS), equation 13.1.1), 



zW" + {b — z)W' — aW = 0, with h = 2 and a = 2a /To- The solutions to Kummer's equation are 
confluent hyper geometric functions, M(a,b,z) and U{a,b,z). The condition that zero torque is 
applied to the disk at the origin now becomes z'^W — > as 2; — > 0. However, as \z\ ^ 0, [/ ~ kl^"'' 
and thus the contribution of the Kummer U function provides a finite torque at the origin. For 
physical reasons we exclude the U function from the solution. Thus the spatial part of W is just 
given by the Kummer M functions. Using the asymptotic expression for M(a,b,z) (AS 13.5.1), 
we find that imposing the condition VF(2:) ^ as z ^ 00 requires that the real part of a satisfy 

< Re(a) < 2 , that is, < Re(cj) < F^ . (10) 
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This constraint on Re(a) is also required if the intercepted flux dL (see equation (2.13) of Pringle 



(1996)1) is to remain finite as z ^ oo. It is not possible to impose the stronger regularity condition 



that \z'^W\ oc \RW\ = \PR\ {i.e., the height of the warp above the untilted equatorial plane) 
remain finite as z — > oo, as z'^W diverges for all Re(a). Real astrophysical disks will always have 
finite RW, however, either due to physical truncation of the disk or because the surface density 
declines to the point where the disk becomes optically thin to the bulk of the incident flux. 

A particular case of interest occurs when Re(a) = 1, that is, a = 1 — iiv, with uj real. This 
corresponds to a steady precession of the disk structure with angular velocity ro/2, and an 
exponential growth of disk structure at a rate e"^*. In this case the M functions are simply related 
to the L = Coulomb functions: 

M(l - iio, 2, 2ix) = e'^'Foiuj, x)x-^/Co{io) (11) 

(AS 13.6.8) where uj can be either positive or negative and Fq is the regular Coulomb wave 
function with L = (AS 14.1.4; note that both uj and x are real). The azimuth of the line of 
nodes is given by 7 = x, i.e., 7 oc i?^/^, which corresponds to a prograde spiral (c/. P96). 

If a; = 0, there is no growth and the disk structure is steady in the precessing frame. In this 
case M simplifies to the solution 

M(l,2,2ix) =e^^^^^ . (12) 

X 

It is interesting to compare this solution (zero torque at the origin) with the torqued solution, 
equation (8), found above. Both solutions precess at the same rate, ro/2, and both solutions 
display an envelope to \W\ of the functional form R^^^^. In both cases, the precession is driven 
solely by radiation (and thus the precession rate depends only on To). However, the zero torque 
solution maintains its constant shape (in the precessing frame) by using viscous forces, which 
would tend to smooth out the oscillations of /3 with radius, to balance the radiation induced 
growth which comes about because d^/dR > (c/. P96, equation [3.10]). 

Some general points regarding the M functions are worth mentioning. As 2; ^ M{a, b,z) ^ 1 
for all a, 6, indicating that the disk has a finite tilt at the origin. Note that since we have required 
Re((T) > 0, the tilt is not stationary in the inertial frame (a = is discussed in section 2.1). As 
we will see below, however, the tilt at R = for the growing modes can be extremely small. 
The Re(a) = 1 M functions are the only genuinely oscillatory (i.e., zero-crossing) functions; the 
radial part \M\ is symmetric about a = 1 for < Re(a) < 2 but the angular dependence differs 
for Re(a) > 1 and Re(a) < 1. Figure 1 shows \M\ oc |/?| as a function of \z\ for several values 
of Re(a) with Im(a) = —uj = 0. For u; > (corresponding to growing modes), the first zero (or 
minimum) of |M| moves to larger \z\, and there is a dramatic increase in the amplitude of the first 
maximum of \M\, which also shifts to larger \z\. Both the magnitude of the shift and the increase 
in amplitude increase with uj. For example. Figure 2 plots \M\ for Re(a) = 1 and several values 
of uj; the behavior for other values of Re(a) is very similar. All the \M\ have been normalized to 
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a maximum of one, and are labeled with uj and the normalization constant /3o- For uj = 10, the 
normalization constant /?o ~ 2.5 x 10^^^. 

As we show below, the maximum value Lt^max of u) (and therefore of the growth rate a) is 
constrained by the outer radius of the disk (either the physical edge or the optically thick edge) . 
There are, however, likely to be unstable modes present throughout the range < u; < Wmax- This 
can be seen by considering a WKB-type analysis and replacing d/dR by ik [cf. P96). In this case 
we obtain 

a = il-TokR+^^ \ . (13) 

Thus for a fixed kR {i.e., a fixed radial mode structure) we see that u oc —\m.{a) is largest, Wmax) 
at large radius. Conversely, as R decreases it is evident that all unstable values of uj in the range 
< u; < Wmax are present. Physically, this arises because the viscosity {v oc R^/'^) has a stronger 
radial dependence than the radiation torque parameter (F oc R). Thus for arbitrarily small oj 
there is always some radius at which the torque term dominates over the viscous term, and the 
disk is unstable to warping. In a real system, however, the development of small-scale modes will 
be limited by the disk age, as we discuss below. 



3. Time Evolution and Steady State 

I 1 

Pringle (1996) concluded that centrally- illuminated accretion disks are unstable to warping 
beyond a critical radius i?c, given by 



Rc ^ / 2^/2^r? y 
Rs \ € ) 



(14) 



the disks remain flat interior to this radius. At first glance this appears substantially different from 
the results of section which showed that the tilt {3 must be non-zero at the origin. However, 
the time evolution of the tilt is crucial. With increasing w, the initial growth of the warp occurs 
at larger radius, and the amplitude of the maximum of [3 relative to the = value grows 
dramatically (c/. Figure 2); for Re(a) = 1, /3max//3o oc (exp(27ru;) — l)/a;. Thus the disk will begin 
to warp initially at the largest possible radius (close to the disk physical or optical depth edge) 
and will be flat interior to this radius, as in Pringle's analysis. The main difference is that here 
the warp has non-zero precession frequency {i.e., the real part of a is non-zero), whereas in P96 
the growing modes were purely imaginary. Note also, from equation (6) for z, that if the warp 
begins at some value of I 

-^1 — I 1 10 5 this corresponds to a. radius 



Rw 1 /?? 



2 



z\l (15) 



Rs 2\e, 

which shows the same scaling with i] and e as equation (4.3) of P96, given above. 
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If we require that the disk be unwarped beyond some maximum radius i?max (as will be the 
case, for example, if the disk goes optically thin beyond some radius, so that the instability ceases 
to operate), then we must have Re(a) = 1, as this is the only mode that passes through zero. The 
warp then has a unique precession rate, Re((7) = Tq/^. Furthermore, the requirement that the 
first zero of (3 (and hence |M|) occur at .Rmax fixes a unique value for the the maximum value of uj 
{cf. Figure 2). (A similar boundary condition will determine Wmax even if a 7^ 1, since the radius 
of the warp onset always increases with oj.) 

One might expect that the steady state of such an accretion disk will be the lo = solution. 
However, real accretion disks will not necessarily reach this state. The tilt /? - and therefore the 
amplitude of the warp - will saturate at some finite value, due to nonlinear effects which have 
been ignored in our analysis. This will of course occur first for the fastest growing mode, at the 
largest unstable radius. With time the slower-growing modes (only those with zeros at i?max if 
we require the disk to be flat there) will develop at smaller radius. However, due to the rapid 
reduction in amplitude /3max//?o with decreasing uj, the slower-growing modes have to evolve 
through a progressively larger number of e— folding times to reach the amplitude of the mode with 

= "^max- For Re(a) = 1, the number of required e— folding times is approximately vr(LiJmax — ^) 
(assuming 2ttuj ^ 1). Thus the minimum value of lo that can have developed significant amplitude 
is set approximately by the condition 

7r(a;max — i^min) / % 
Tdisk = Ti(^maxj (loj 

where Tdisk is the lifetime of the accretion disk and rj(a;max) is the timescale for growth of the 
mode with m = Wmax- As we will see in section this is likely to be pose a significant constraint 
on Wmin for accretion disks around supermassive black holes in active galactic nuclei. 

In addition, the effects of shadowing - which we have ignored - may be important for modes 
which have multiple zeros between R = and i?max- This will depend on both the actual size of 
the central source of luminosity and the amplitude at which the warp eventually saturates. 

Figure 3 shows the shape of the warp for the a = 1, w = (steady-state) solution. The disk 
surface has been plotted out to the radius (x = 2tt) at which the warp returns to the equatorial 
plane. The maximum amplitude of the warp has been fixed at 10% of this radius. The spiral 
shape of the line of nodes is readily visible. 



4. Applications 

4.1. The Warped Maser Disk of NGC 4258 

VLBI observations of the 22 GHz water maser line emission from the nucleus of NGC 4258 
(which possesses a low-luminosity active nucleus) have shown that the maser emission arises in a 
thin, warped disk, with an inner radius of 0.13 pc and an outer radius of 0.25 pc, which exhibits 
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a perfectly Keplerian velocity curve to within the measurement errors ( IGreenhill et al. (1995)| 



Miyoshi et al. (1995)). The derived central mass is = 3.6 x lO'^ -A^©, and as shown by [Maoz 
1995), dynamical arguments essentially rule out mass models other than a supermassive black 



hole. The obliquity of the warp is ~ 0.25 at the outer edge of the masing zone (Greenhill 1995, 
private communication) and decreases steadily inward. 

I , I 

Neufeld, Maloney, & Conger (1995) showed that X-ray irradiation of dense molecular gas is 

an effective mechanism for generating powerful water maser emission; this also naturally explains 
the association of water megamasers with active galactic nuclei ( Claussen &: Lo (1986)| ; BraatzJ 
Wilson, fc Henkel (1994)| ). These models also show that there is a critical pressure at which 
the gas goes from molecular to atomic. By modeling the NGC 4258 disk using the standard 
a— prescription for viscosity, Neufeld &: Maloney (1995)| (NM95) showed that the disk would 
be extremely thin {h/R <^ 10~^), that the transition from molecular to atomic gas should be 
identified with the outer edge of the masing region, and that the mass accretion rate through the 
disk is M/a ~ 7 x 10^^ Mq yr~^. This mass accretion rate in turn implies that the radiative 
efficiency e ~ 0.1. The existence of the warp in this picture is crucial, since for a flat disk 
neither the greatly reduced irradiation from the central source nor heating by viscous dissipation 
would keep the disk warm enough to produce substantial amounts of water or excite the masing 
transitions. In fact, NM95 suggested that the inner edge of the masing disk marked the inner 



boundary of the warp, as is the case in the model for the warp derived by Miyoshi et al. (1995) 



Although an arbitrary imposed warp would be kinematically stable in a Keplerian potential, since 
the vertical and azimuthal frequencies are the same, the origin of the warp in NGC 4258 remained 
puzzling. 



Here we show that the radiation-driven warping instability discovered by [Pringle (1996) 



provides a plausible explanation for the warping of the disk. From equation (6) for z, with 
Rg ~ 1.1 X 10^^ cm for the central object in NGC 4258, we find that the inner and outer edges of 
the masing region are at 

l^lin « 28eo.i , l^lout ~ 39eo.i (17) 

where e = O.leo.i- As can be seen from Figure 2, this requires Wmax ~ 10 for eo.i ~ 1, as is indicated 
by the results of NM95. (We note that in the molecular zone at the outer edge of the masing 
region in the model of NM95, the hydrogen column density through the disk Nh ~ lO'^^ cm~^, 
satisfying the assumption that the disk is optically thick to both absorption and emission. Also, 
the outer edge of the masing zone does not have to coincide with the maximum of the warp.) It is 
likely that Wmin is not much smaller than this value, due to the limit set by the finite age of the 
accretion disk. The characteristic timescale for growth of the instability is 

niu) ^1- = ^ (18) 

\a\ iOL o 

assuming that to ^ Re(a). This can be expressed in terms of the viscous timescale 
Tvisc ~ 2'kR^T,/M: 

12 

Ti{i0) r |2:r-^r?"Vvisc (19) 

U) 
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which, under the assumption that the luminosity is accretion-powered, depends only on e and 
not on L and M individually. Assuming standard a— viscosity and negligible self-gravity for the 
accretion disk, this becomes 

, , 4 cRs 7 Ma , , 

n{uj) f « 3 X 10' yrs (20) 

oj aCgC cjai 360.1 

where Mc = 10^ Mg Mq, Cg is the sound speed, the gas temperature T = lO^Ts K, and the 
numerical coefficient is for molecular gas. In the masing molecular zone, 0.3 < T3 < 1. While 
the characteristic growth time for the a; = 10 mode in NGC 4258 is Tj ~ 5 x 10^ years, the time 
required for even the w = 5 mode to become comparable in amplitude is an order of magnitude 
longer. Thus we expect that the warp will be dominated by modes with oj close to Wmax, unless 
the accretion disk is extremely long-lived (rjisk 3> 10^ years). Hence the predictions of the 
radiation-driven warp instability model are in reasonable agreement with the observations of NGC 
4258, for e ~ 0.1. 

This instability cannot provide an explanation for the warp if NGC 4258 possesses an 
advection-dominated accretion disk with e ~ 10~^, as suggested by Lasota et al. (1996) . From 



equation (17), this requires that Izjin, l^^lout ^ 1- However, the only modes that have any power 
at such small \z\ have to <\. Equation (20) shows that the timescale for growth of these modes is 
prohibitively long, with Tj > 10^ years. 



4.2. The 164-day Precession Period of SS 433 

SS 433 is an eclipsing binary system containing an early-type star and an accretion disk around 
a compact object. The optical and infrared spectrum shows two sets of highly Doppler-shifted 
emission lines, which reach velocity extremes of -|- 50, 000 and —30,000 km s"^, respectively. These 
emission lines are believed to arise from oppositely-directed relativistic {v/c ~ 0.26) beams of 
material emanating from SS 433 (reviews of observations and models of SS 433 can be found in 



Margon (1984) and Zwitter et al. (1989) ). The high- velocity systems show systematic changes in 
wavelength, with a period of approximately 164 days. This period is generally regarded to be the 
result of precession of the ejection axis of the system. There is also direct evidence for precession 



of the beams from radio maps ( [Hjellming &: Johnston (1986) ). In addition, there is evidence for 



precession of the accretion disk, which dominates the optical light of the system, from optical 
photometry ( |Margon (1984) ). 



A number of models have been proposed to explain the precession of the accretion disk. Most 
of these involve "slaved" precession, in which the residence time in the disk is sufficiently short 
that the inner disk follows the precessional motion of a misaligned companion star or outer disk 
{e.g., Katz et al. 1982). This in turn requires that the disk be very thick, if the internal disk 



pressure is to be comparable to the very large viscous stresses which are required. Support for the 
existence of a thick {h/R ~ 1) disk is provided by analysis of the light curve ( Anderson, MargonJ 
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& Grandi (1983)). 

However, if the accretion disk in SS 433 is substantially warped, the actual thickness of the 
disk may be much smaller than inferred from the light curve, which will measure the amplitude of 
the warp. (This is certainly the case in NGC 4258, in which the solid angle occupied by the disk 
about the source is completely dominated by the warp.) As shown in section |2|, the real part of 
a is in general not zero, so the warp will precess with a fixed pattern speed. Assuming that the 
accretion disk in SS 433 is subject to the radiation-driven warp instability, we now show that that 
the 164-day precession period can be explained as the result of precession of the warp. 

We again assume that Re(a) = 1, so that we can require the outer boundary of the disk to be 
unwarped. As noted earlier, this fixes the precession rate at To/2 and so the precession timescale 
is just Tprec ~ 'i'^/To- Since the assumption of an isothermal disk is likely to be poor for SS 433, 
we instead write the viscosity in the form 



a 



■" = 7i^ni:j 111 <^'> 



(22) 



and so the precession period is 

ae c \h J 

where this should be evaluated at R = i?warp- For accretion disks around Mc ~ solar-mass objects, 
the timescales for growth of the unstable modes are sufficiently short that the warp reaches the 
steady-state, uj = solution. For a = 1 the first maximum occurs at \z\ ^ 1.3tt, and we take this 
to correspond to i?warp- Thus the precession timescale becomes 

rprec~ 430^— - . 23 

For SS 433, e is probably <C 0.1, since the accretion rate is almost certainly substantially 
super-Eddington {e.g., Begelman et al. (1980)| ) and so the precession timescale is 



rf ( Mc \ (R 



which can explain the 164-day precession period in SS 433 for reasonable values of a and {R/h). 
We also note that the precession can be either prograde or retrograde with respect to the direction 
of disk rotation, as can be seen from the derivation leading to equation (2.18) of P96. 



5. Summary 



Our global analysis confirms the conclusion of Pringle (1996) that thin, centrally- irradiated 
accretion disks are generically unstable to warping driven by radiation pressure. The maximum 
growth rate is set by the largest unstable radius (either the physical boundary of the disk or the 
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radius at which it becomes optically thin to the bulk of the incident or re-emitted flux). The 
growth rate on smaller scales in accretion disks around supermassive black holes in active galactic 
nuclei will in general be limited by the lifetime of the accretion disk. Although the final state of 
the instability will be determined by nonlinear processes which have not been considered in this 
analysis, we have shown that the warping instability provides an explanation for the existence 
of the warp in the masing accretion disk in NGC 4258. Furthermore, for plausible parameters 
the 164-day precession period in SS 433 can be identified with the precession period of the warp, 
as can the long (30 — 300 day) periods seen in other X-ray binary systems, such as Her X-1 (see 
White, Nagase, &: Parmar 1995| for a review). 

Radiation-driven warping and precession thus offers a robust mechanism for producing 
warped, precessing accretion disks in AGN and in accreting binary systems such as Her X-1 and 
SS 433. Because the warp adjusts itself in such a way as to precess at a uniform rate at all 
radii, radiation-driven warping is an inherently global mechanism, thereby avoiding the difficulties 
inherent in other proposed mechanisms for producing warping and precession. This instability 
should generally be important in thin accretion disks. 

PRM acknowledges support through the NASA Long Term Astrophysics Program under 
grant NAGW-4454. MCB was supported by NSF grant AST 91-20599. 
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Fig. 1. — The absolute magnitude of the Kummer function M(a, 2, z) for u = — Im(a) = and 
several values of Re(a) (i.e., the time- independent solutions of the twist equation). The magnitude 
of the local angle of tilt \f3\ is equal to \M\ times a normalization coefficient. The radial variable 
z is defined hy z = iV2e{R/Rs)^^'^ /r], where e is the radiative efficiency e = L/Mc^, Rs is the 
Schwarzschild radius, and r) is the ratio of vertical to radial viscosities. The Kummer function with 
Re(a) = 1 is the only zero-crossing function. For other values of a \M\ is symmetric about a = 1. 
Requiring that the tilt /? — as i? — oo requires that the real part of a satisfy < Re(a) < 2. 

Fig. 2. — The amplitude of the local angle of tilt for several of the growing solutions of the twist 
equation. The Kummer functions |M(1 — iuj,2,z)\ are plotted. All these modes have Re(a) = 1 
and grow as e'^*. They have all been normalized so that the maximum amplitude is 1; each curve is 
labeled with to and with the normalization constant (3o, i-e., the magnitude of the tilt at the origin. 
For a; ^ 1 this is extremely small, so that for modes with high growth rates the disks remain flat 
interior to a critical radius. The maximum value of to for a disk is fixed by the largest unstable 
radius, which is either the physical edge of the disk or the radius at which it becomes optically 
thin. 

Fig. 3. — Surface plot of the warp shape for the steady-state (a; = 0), a = 1 mode of an isothermal 
accretion disk. The warp has been plotted out to the radius at which it returns to the equatorial 
plane; the amplitude of the warp has been fixed at 10% of this radius. 
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